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1. Introduction 

A remarkable feature of the v\' meson in comparison to the pions, the kaons and the r\ meson is 
its extraordinary high mass of about 1 GeV: while pions consist only of light quarks and - therefore 
- exhibit a rather small mass around 140 MeV, strange quark contributions rise the mass of the four 
kaons and the r\ meson to 500 - 600 MeV. However, the larger strange quark mass cannot explain 
the mass value of the r\' meson. On the QCD level, the reason for the mass gap is thought to be the 
breaking of the Ua ( 1 ) symmetry by quantum effects. 

Although lattice QCD provides access to flavour singlet pseudo-scalar states, significant con- 
tributions from quark disconnected diagrams complicate the determination of their properties. This 
might be the reason for the rather short list of publications covering these mesons. Recent studies 
for 2 + 1 dynamical quark flavours can be found in [jj], ^, |3]] and pj|. For a recent publication with 
Nf = 2 Wilson twisted mass fermions see Ref. [Q]. 

In this proceeding contribution, we build on a recent paper of some of the authors ^ and 
discuss the determination of 17 and T]' meson masses and the mixing angle using the Wilson twisted 
mass lattice QCD formulation with Nf = 2 + 1 + 1 dynamical quark flavours. Compared to 
Ref. [^] we present results for an additional ensemble, which allows us to better control systematic 
uncertainties. The final results stay virtually unaffected compared to Ref. 

The results presented in Ref. [^] were extracted for the unitary case, with identical sea and 
valence quark regularisations. In addition we discuss a mixed action approach with so called 
Osterwalder-Seiler strange and charm valence quarks [||]. Possible advantages of such an approach 
are a powerful variance reduction technique for strange and charm quarks Q5| ^] and reduced twisted 
mass induced isospin splitting in the strange/charm doublet. We discuss how to perform the match- 
ing of valence and sea formulation and present first results for the variance reduction. 



2. Lattice actions 



The results we present are based on gauge configurations provided by the European Twisted 
Mass Collaboration (ETMC) and correspond to three values of the lattice spacing, a = 0.061 fm, 
a = 0.078 fm and a = 0.086 fm. The pion masses range from 230 to 500 MeV [M O]. A list of the 



investigated ensembles is given in Table [T[ For setting the scale we use throughout this proceeding 
contribution the Sommer parameter ro = 0.45(2) fm [1 1 ]. 
The Dirac operators for u and d quarks [^] reads 



Dt=D w + m + iH(Y 5 i; i 



(2.1) 



with \i( the bare light twisted mass parameter. For the heavy doublet of c and s quarks [ |12| ] the 
Dirac operator reads 



D h =D w + m + /jU^t 1 + M^ 3 



(2.2) 



Dyy denotes the standard Wilson operator and the T ! are Pauli matrices acting in flavour space. 
The value of the bare quark mass mo was tuned to its critical value [13, [H]], leading to automatic 
order G (a) improvement at maximal twist [14], which represents the most notable advantage of 
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tmLQCD. The two bare quark masses }X a , n§ are related to the physical charm and strange quark 
mass via 

m c>s = n a ± Zn s , (2.3) 

where Z = Zp/Z$ defines the ratio of pseudo-scalar and scalar renormalisation constants Zp and 
Z5. We denote quark doublets in the physical basis by Ye,h and in the twisted basis by Xe.h- In the 
continuum they are related by exact axial rotations 

W,h = e in ^ 4 X e,h , %,h = Xl,he inrst3 ' l/4 ■ (2-4) 
The main drawback of this formulation is the breaking of flavour symmetry at finite values of 



the lattice spacing, which was shown to affect mainly the value of the neutral pion mass [15, 16, 



170. Furthermore, for the non-degenerate quark doublet this introduces mixing between charm and 
strange quarks. 

This complication can be circumvented in a mixed action approach using so-called Osterwalder- 
Seiler (OS) type quarks in the valence sector for strange and charm quarks [§J. Formally, we in- 
troduce one twisted doublet each for valence strange and valence charm quarks [|| |^]. Therefore, 
in the valence strange and charm sector flavour mixing is avoided. It was shown in Ref. []|] that 
automatic G{a) improvement is not spoiled by this approach and that unitarity is restored in the 
continuum limit. 

Formally, we have introduced two strange (charm) quarks differing in the sign of the twisted 
mass. We will denote the one with positive sign with s (c) and the other one with s' (c'). With jj. s 
and ;U C we denote the bare OS strange and charm twisted masses, respectively. 

The two actions in the sea and the valence sector must be matched appropriately. This match- 
ing can be performed using different observables: firstly, one can match using the quark mass 
values 

Ms,c = Ma ± Z fig , 

requiring the knowledge of the ratio of renormalisation constants Z; secondly, unitary and OS kaon 
mass values can be matched. For the OS kaons we have the possibility to define the following OS 
interpolating fields in the twisted basis 

0&* = Xs> Id (x) , K+ = %s ijsXd , 

which will lead to different kaon mass values at finite values of the lattice spacing. The correspond- 
ing kaon masses will be denoted by M K os and M K + , respectively. The unitary kaon mass value is 
denoted by Note that there is no isospin splitting for the kaons in Nf = 2 + 1 + 1 Wilson 



twisted mass lattice QCD [19]. 



Thirdly, one can use the mass of the artificial rj s meson - a pion made out of strange quarks 
- to match sea and valence actions. can be determined from the connected only correlation 
function of the following operator 

= -^.XsXs-Xs'Xs')- 

Similarly, one can define matching observables for the charm quark. However, as we are 
mainly interested in 17 and tj' mesons, we expect little impact from the charm quark. Therefore, we 
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ensemble 


P 


a}X( 


a\i a 




L/a 


Aconf 


N s 


N b 


A30.32 


1.90 


0.0030 


0.150 


0.190 


32 


1367 


24 


5 


A40.24 


1.90 


0.0040 


0.150 


0.190 


24 


2630 


32 


10 


A40.32 


1.90 


0.0040 


0.150 


0.190 


32 


863 


24 


4 


A60.24 


1.90 


0.0060 


0.150 


0.190 


24 


1251 


32 


5 


A80.24 


1.90 


0.0080 


0.150 


0.190 


24 


2449 


32 


10 


A 100.24 


1.90 


0.0100 


0.150 


0.190 


24 


2493 


32 


10 


A80.24s 


1.90 


0.0080 


0.150 


0.197 


24 


2517 


32 


10 


A 100.24^ 


1.90 


0.0100 


0.150 


0.197 


24 


2312 


32 


10 


525.32 


1.95 


0.0025 


0.135 


0.170 


32 


1484 


24 


5 


535.32 


1.95 


0.0035 


0.135 


0.170 


32 


1251 


24 


5 


555.32 


1.95 


0.0055 


0.135 


0.170 


32 


1545 


24 


5 


575.32 


1.95 


0.0075 


0.135 


0.170 


32 


922 


24 


4 


585.24 


1.95 


0.0085 


0.135 


0.170 


24 


573 


32 


2 


D15.48 


2.10 


0.0015 


0.120 


0.1385 


48 


1045 


24 


10 


D30.48 


2.10 


0.0030 


0.120 


0.1385 


48 


469 


24 


3 


D45.32sc 


2.10 


0.0045 


0.0937 


0.1077 


32 


1887 


24 


10 



Table 1: The ensembles used in this investigation. For the labeling we employ the notation of ref. [|To|]. Ad- 
ditionally, we give the number of configurations N con f, the number of stochastic samples N s for all ensembles 
and the bootstrap block length Nt,. The D30.48 ensemble was not yet included in Ref. []|]. 



will not investigate different matching observables for the charm quark and use only one /^.-value. 
It is obtained by matching aM K + = aM^ to determine a^L s . The latter is then used to obtain ajx c via 
the relation Eq. (2.3). 

For the extraction of quark masses and pseudo-scalar decay constants it turned out that using 
M K + for the matching lead to the smallest lattice artifacts [|20|]. However, this may depend on the 
quantities under consideration. 



3. Pseudoscalar flavour-singlet mesons 

In order to study properties of pseudoscalar flavour-singlet mesons we have to consider light 
strange and charm contributions to build a suitable correlation function matrix. Hence, in the 
physical basis we are after computing 



( m(t)i\ t (o) m{t)r] s (o) t]KO^(o) n 

r\s(t)nt(0) Vs(t)riM i?,(/)rj c (0) | , 
\r]c(t)r]m r} c (t)r} s (0) J] c Wt? c (0), 



involving the following interpolating operators 

m = (Yui75Yu + W75Wd)/y/2, rj s = (\j/ s ij 5 Y s ), r\ c = (WftVc)- (3-2) 
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Since we are working in the twisted basis we have to rotate these operators. Let us start with the 
unitary action, where the rotation in the light sector is given by [||] 



^2 ' XuXu ~\~ Xd Xd ) 



(3.3) 



where the left- and right-hand side correspond to physical and twisted basis, respectively. Also for 
strange and charm we have to consider a doublet as follows 



iY5- 



1±T J 




(3.4) 



The flavour space projector ( 1 ± T 3 ) /2 distinguishes between charm and strange contributions in the 
physical basis. In the twisted basis we need to consider the following operators when calculating 
correlation functions 



Z(xd75Xc - XsijsXs)/ 2 ~ (XsXc + XcXs)/2 , 
Z{xJY5Xs ~ XdY5Xc)/2 - {XsXc + XcXs)/2 ■ 



(3.5) 



Note that the ratio of renormalisation constants Z appears in the sum of pseudoscalar and scalar 
currents. However, Z is not needed for extracting the masses of f] and r\' , as we explain in Ref. [§]. 
For further details on how to construct the correlation matrix we also refer to Ref. [§]. 

Let us now discuss the rotation for the mixed action approach with OS valence strange and 
charm quarks. For the light quarks the operator is identical to the unitary one in Eq. (3.3), i.e. 
<^? s = @g. But for strange and charm quarks the operators are significantly simpler, due to no 
flavour mixing in between strange and charm: 



1 . _ . , 1 



Xc'Xc 



c 



(3.6) 



Again, we build a correlation function matrix of the form given in Eq. (3.1). Obviously, there is no 
mixing of pseudoscalar and scalar currents like in Eq. (3.5) and, therefore, no ratio of renormalisa- 
tion constants appears. 

For both, the unitary and the mixed action approach we solve the generalised eigenvalue prob- 
lem nng 

V{t) nM(t,t ) = Xto(t 1 to) n*o) T7 (,!) (Mo)- (3-7) 
Taking into account the periodic boundary conditions for a meson and solving 

A(»)(Mo) e -m^t + e -m(")(T-t) 



(3.8) 



AW(/ + l,?o) e -m<»)(?+l) +e -m(»)(r-(f+l)) 

we determine the effective masses m^"' , where n counts the eigenvalues. The state with the lowest 
mass should correspond to the 77 and the second state to the 77' meson. Alternatively, we use a 
factorising fit of the form 



2mW 



exp(-m w ?) + exp(-m w (T-t)) 



Mi 



(3-9) 
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to the correlation matrix matrix % '. The amplitudes A qjl correspond to (0\qq\n) with n = 77 , 77', ... 
and q = £,s,c. Note that for this physical interpretation of the amplitudes the ratio of renormalisa- 
tion constants Z is unavoidably required as input for the unitary approach 

3.1 Variance reduction 

In general the correlation functions consist of quark connected and disconnected diagrams. 
The connected pieces have been calculated via the so called "one-end-trick" [ p4| ] using stochastic 
timeslice sources. For the disconnected diagrams we resort to stochastic volume sources with 
complex Gaussian noise. The light disconnected contributions can be estimated very efficiently 
using the identity ^ 

D- 1 - D- 1 = -2wD- 1 75 D- 1 . (3. 10) 

In the heavy sector of the unitary setup such a simple identity does not exist. Instead we use the 
(less efficient) so called hopping parameter variance reduction, which relies on the same identity 



as in the mass degenerate two flavour case (see ref. [24] and references therein). The number of 
stochastic volume sources N s per gauge configuration we used for both the heavy and the light 
sector is given for each ensemble in table |l[ In order to check that the stochastic noise introduced 
by our method is smaller than the gauge noise we have increased N s from 24 to 64 for ensemble 
#25.32, which did not reduce the error on the extracted masses. 



For OS strange and charm quarks the variance reduction trick Eq. ( 3.10| ) is also applicable, as 



noted in Ref. [gj. Consider the disconnected contributions for the operator in Eq. (p.q). Again, we 
may write for the strange quark 

D- 1 - D-, 1 = -2i/x s D;, 1 75 D~ 1 , (3.11) 

and similarly for the charm quark. This identity can be used like in the light sector to compute the 
disconnected contributions of strange and charm quarks to the OS correlator matrix with greatly 
reduced noise [HI. 



4. Results from Unitary Strange and Charm Quarks 

Most of the unitary results have already been published in Ref. [§], the only exception is 
ensemble D30.48. However this additional point does only very mildly affect the final results. 

We have calculated all required contractions for the correlator matrix using local and fuzzed 
operators, yielding a 6 x 6-matrix. The number of gauge configurations per ensemble is given in 
Table |. All errors have been calculated from bootstrapping with 1000 samples. To compensate 
for autcorrelation we have used blocking, the number of configuration per block Nb is also given 
in Table [l] and it was chosen such that the resulting blocklength in HMC trajectories A^hmc fulfils 
Nhmc > 20. For a more detailed discussion on autocorrelation, which mainly affects the T] '-state 
we refer to [Q]. 

4.1 Extraction of Masses 

The details of our GEVP and fitting procedures to extract f] and f]' masses are explained in 
Ref. [|6p. In Figure [T] we show the masses of the r] (filled symbols) and Tj' (open symbols) mesons 
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A-Enscrnblcs 
A80.24s, A100.24s 
B-Enscrnbles 
D-Ensomblcs 
D45.32sc 



0.2 0.4 0.6 0.8 1 1.2 1.4 

\2 



(a) 



1.5 




0.5 



constant fits 
linear fits 
A40.32 
A60.24 
B35.32 
B55.32 
D30.48 
D45.32sc 



-B- 



0.01 0.02 0.03 0.04 0.05 

(b) 



Figure 1: (a) 77 (filled symbols) and 77' (open symbols) masses in units of ro- (b) roM^ as a function of 
(a/ro) 2 for the ensemble sets Si and S2. 



for the various ensembles we used as a function of the squared pion mass, everything in units of 
ro. The values of the chirally extrapolated for each value of j3 can also be found in Ref. [|]]. 
We present the values for aM^ and aM^i together with kaon and pion mass values in Table ^| with 
statistical errors only for the two ensembles D30.48 and D45.32sc. The results for D30.48 are new 
compared to what was shown in Ref. [Q]. D45.32sc will be used for the mixed action analysis. All 
other results can be found in Ref. [§]. 

It is clear from the figure that the 77 meson mass can be extracted with high precision, while 
the 77' meson mass is more noisy. 



4.2 Scaling Artifacts and Strange Quark Mass Dependence of M n 

The results displayed in the left panel of Figure [j] have been obtained using the bare values of 
afi a and apL§ as used for the production of the ensembles. Those values, however, did not lead to 
the physical values of, e.g., the kaon and D-meson masses JIo| , |I9[ ]. 

For M r] the statistical uncertainty is sufficiently small to attempt to correct for the mismatch 
in the strange quark mass value and to try a scaling test. For this we need to compare at the 
three different values of the lattice spacing for fixed values of for instance toMk, toMd, roMps and 
the physical volume. From volume and the charm quark mass value we expect only little influence 
given our uncertainties and hence, we are going to disregard these minor effects in the following. 

As discussed in Ref. [g], we have to perform an interpolation of in Mr. For this purpose, 
we treat the masses of the 77 -meson and the kaon like in chiral perturbation theory as functions 
M 2 =M 2 [Mp S ,Mj|] and define the dimensionless derivative 



'd(aMn) 21 



(4.1) 
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1.5 



0.5 




^ 0-^^77,exp 

A-Enscmblcs 
A80.24s, A100.24s 
B-Enscmbles 
D-Enscmblcs 
D45.32sc 



0.2 0.4 0.6 0.8 1 

\2 



1.2 1.4 



(r M PS ) z 



(a) 



1.5 




0.5 



(MJMk)^ >— ■- 

{AIJAI K )nh ys ' B- 

A-Ensembles 1 — ^- 
A80.24s, A100.24s 1 — ■- 
B-Ensembles 1 — 
D-Ensembles 1 — t 
D45.32sc 1 — ♦- 



0.2 0.4 0.6 0.8 

(r () A/ps) 2 

(b) 



1 



1.2 1.4 



Figure 2: (a) r(,M n as a function of (r M PS ) 2 . (b) M^/M K as a function of (r M PS ) 2 



Next we make the approximation that is independent of the quark mass values iXi,ii. a ,jXg and 
j8. Its value we can estimate from A80.24 and A%Q.2As as well as from A100.24 and A100.24s. On 
average we obtain D r] = 1.60(18). 

Now we use this value of D r] to correct two sets S\,S2 of three ensembles, namely Si = 
{A40.32,B35.32,D30.48} and S 2 = {A60.24,B55.24,D45.24} to a common value of r M K « 1.34 



using 



(r M TJ ) 2 = (r M IJ ) 2 +D TJ -A K , 



where Ak is the difference in the squared kaon mass values to the squared reference values (in units 
of ro). For each set the three points have approximately fixed values of roMps. 

We plot the resulting r^M-q values for sets 5^2 as a function of (a/ro) 2 in the right panel 
of Figure fil. Both data sets are still compatible with a constant continuum extrapolation giving 
rnA ^Si .const = 1-447(45) and r Q M^ const = 1.480(34), respectively, which we indicate by the 
horizontal lines. We can also perform a linear extrapolation, leading to roM"^ lin = 1.60(25) 
and ro^^, ii n = 1-61(14), which is also shown in the figure. The difference in between the two 
extrapolated values for each set are 



r AM! 



J?, Si 



0.15(25), r AM: 



a^0 



0.13(13) 



(4.2) 



and they give us an estimate on the systematic uncertainty to be expected from the continuum 
extrapolation. Both results agree well, although the one for Si exhibits twice the error. We will 
therefore quote an 8% relative error from AM^ 2 /M^° const for our mass estimates, which was 
already used in Ref. [|]] where S\ was not yet available. 

In order to obtain a more complete picture, we now correct all our ensembles for the slightly 
mistuned value of Mr. For this we follow the procedure which was discussed in detail in [Sj, i.e. 
we shift the kaon mass values for all ensembles to a common line (roM/f) 2 [(roMps) 2 ] determined 
such that it reproduces the physical kaon mass value at the physical point. Next we correct the rj 
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ensemble 


aMps 


aM K 




aMrji 


D30.48 


0.09776(45) 


0.17760(23) 


0.205(16) 


0.38(4) 


D45.32sc 


0.07981(30) 


0.17570(84) 


0.192(15) 


0.30(4) 



Table 2: Results of aM-q, aM^i for ensembles D30.48 and D45.32.sc. and the corresponding values for the 
charged pion mass Mps and the kaon mass Mk- The D30.48 is new compared to Ref. [||] and ensemble 
D45.32sc is needed for the mixed action results. Results for all the other ensembles can be found in Ref. [Hi. 



masses appropriately. The result of this procedure is shown in the left panel of Figure En we show 
values of the corrected r\ masses r^M-q for all our ensembles as a function of (roMps) 2 . It is evident 
that all the data fall on a single curve within statistical uncertainties, which confirms that M-q is not 
affected by large cut-off effects. Note, however, that we again ignored possible \in, \i c , jig and /3 
dependence with this procedure. 

4.3 Extrapolation to the Physical Point 

Since we have now fixed the strange quark mass to its physical value using M^f = 498 MeV, 
we can attempt a linear fit to all corrected data points for (roM^ ) 2 [(roMps) ] . Using ro = 0.45 (2) fm 



as in Ref. [JTTJ], the fit yields r M,, [r^M 2 ] = 1.256(54) stat (100) sys and in physical units 

M r] (M x ) = 55 1 (33) stat (44) sys MeV , (4.3) 

where the experimental mass-value of the neutral pion M n o = 135 MeV has been used for M K . In 
the SU (2) chiral limit we obtain r M° = 1.230(65) sta t(98) sys or M° = 539(35) stat (43) sys MeV . 

As the procedure used to correct the r\ mass for mistuning of the strange quark mass ignores 
a possible dependence on pL a , jJ.$ and j3, it is desirable to have a cross-check. In Ref. [^] we 
discussed two possible options for this. The first one is to study the ratio M^/Mk, for which it 
was shown that most of the strange quark dependence cancels. The second possibility is to study 
the Gell-Mann-Okubo (GMO) relation which is motivated by chiral perturbation theory. Here we 
simply repeat the analysis for both cases, including the new data point for D30.48. 

For the ratio (M^/Mk) 2 a linear extrapolation for all available data in (roMps) 2 to the phys- 
ical pion mass point yields (see the shaded band in the right panel of Figure |2|) (M^/Mk) m = 
1.123(26) , which agrees well with the experimental value (A/jj/Mk) = 1.100. Using the exper- 
imental value of M k q we obtain M-q = 559(13) sta t(45) sys MeV. Note that in this analysis the scale 
ro = 0.45(2) fm is only required for determining the physical pion mass point. As the slope of the 
extrapolation is rather small, the statistical uncertainty in is significantly smaller than for the 
direct extrapolation of (roM^) 2 . 

Considering the GMO relation and again performing a linear extrapolation in (roMps) in- 
cluding all available data we obtain (3M 2 / (4Af| - M 2 )) M = 0.970(47) at the physical pion mass, 
which is in agreement with experiment, (3M 2 /(4Mj| — M 2 )) exp = 0.925. Using the experimen- 
tal values of M K o and M K o we now obtain = 561(14) stat (45) sys MeV, where the first error is 
statistical and the second systematic estimated again from the scaling violations discussed above. 
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O.O.N 



0.06 



t 0.04 



0.02 





Figure 3: We show (M„ s ) 2 , {M K o& ) 2 and (M K + ) 2 extracted with OS -type valence quarks as functions of the 



bare OS strange quark mass ajx s . As horizontal lines we show the unitary values of (aMx) 2 and (aM„ mt ) 2 



5. Results using OS strange and charm quarks 

In this section we will discuss first results for v\ and v\' meson masses using OS strange and 
charm quarks. We focus here on one ensemble, namely D45.32sc, see Table @. 



For given values of ajj,( ,ajj, s and a\i c we determine the correlation matrix Eq. (3.1) using the 



interpolating operators Eq. ( |3.3D and Eq. ( |3.6| ) using 900 configurations. Two consecutive config- 
urations are spaced by 4 trajectories of length 1. Disconnected contributions are estimated using 
32 Gaussian volume sources per gauge configuration. We use local and fuzzed operators and build 
correspondingly a3x3or6x6 correlation function matrix. Errors are estimated using 1000 boot- 
strap samples. The data is blocked in blocks of length 5 to account for autocorrelations. Mass 
values of X] and r]' are extracted using the same methods as in the unitary case. 

5.1 Matching The Strange Quark Mass 

As discussed in Section sea- and valence-actions can be matched using various different 
observables. We will consider here M K os, M K + and M^ s . In Figure || we show (aM K os) 2 , (aM K +) 2 
and (aMjjJ 2 as functions of apL s for ensemble D45.32SC. We also show both (aM™ 1 ) 2 and (aM^) 2 
as horizontal lines. The value of aMk can be found in Table ||. The value for aMJj mt has been 
determined using the connected contributions to the unitary correlator matrix only, and its value is 
aM™ x = 0.2105(14). 

From Figure ^] it is first of all clear that using different valence quantities leads to very different 
matching values for afl s . Using for instance M K + leads to a\l s = 0.0149(3), while matching M K os 
leads to a\x s = 0.0102(3). Therefore, for this proceeding we decided to use the following three 
values for the bare OS strange quark mass 

an s = 0.01 , 0.0125 , 0.018 , 0.025 . 

These four values bracket the three matching values stemming from the different matching quan- 
tities. It is worth noting that the leftmost data points in Figure |3| correspond to the case where 
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Figure 4: M° s and M° s as a function of afl s for D45.32.sc. In addition we show as horizontal lines the 
corresponding unitary masses with errors as shaded band. Note that M-q and M^i are based on a larger 
number of configurations. 



ajj, s = ajj,£. Hence, aM K os = aM^. The splitting in between aM K os and aM K + is an €?(a 2 ) effect, 
which disappears in the continuum limit. 

As mentioned before we do not expect the precise charm quark mass value to be important 
for r\ and r]' , which is also confirmed by our unitary results. Therefore, we use only one value 
a\L c = 0.172 for the OS bare charm quark mass. It is obtained by matching <zM K + = aMk to 



determine a\i s . a\i c is then obtained via the relation Eq. ( |2.3| ) 



5.2 The valence strange quark mass dependence of the OS r/ states 

In Figure [| we show the lowest two states extracted from a 6 x 6 matrix as a function of the 
OS strange quark mass a\i s for the DA'S. Use ensemble. The state with the lowest mass should 
correspond to the 77, the second to the r/' state. In addition to the OS results we show as horizontal 
lines the unitary mass values discussed in the previous section (see Table |J). 

In general one observes from Figure £| that the r/ mass value can be extracted with high statisti- 
cal accuracy, while the r\' suffers from similar noise as observed for the unitary r/'. In fact, for the 77 
meson mass we obtain slightly better accuracy as compared to the unitary case. However, the noise 
reduction trick applicable for OS strange and charm quarks does not seem to help for the extraction 
of the OS f]' meson mass: in the effective mass plots a plateau is only hardly visible. Hence, we 
expect that the OS r\' is affected by similarly large systematic uncertainties as the unitary tj'. 

Therefore, we consider in the following only the OS 77 meson and its /i s dependence. It is 
interesting to understand which of the matching observables discussed above yields the best agree- 
ment in between OS and unitary r\ mass values. First of all, we observe a rather mild dependence 
of M° s on the bare OS strange quark mass ajX s , see Figure |] In order to compare to the unitary 
case, we compute 

D ° s W = a8(1) 
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Figure 5: (a) We show the OS to unitary ratio of relative errors of the 77 -eigenvalue as a function of t/a. (b) 
the same as (a) but for the effective masses. 



with statistical error only, and compare to = 1.60(18) defined in Eq. |4.l[ The OS value is 
significantly smaller which we attribute to the large sea quark contributions to the 77 meson mass. 

As a consequence, the OS 77 mass agrees within errors with the unitary 77 mass value for all 
apt s values considered. The agreement is best around the M K + matching point, though. Given the 
large uncertainties in the OS t]' meson mass, we find also for the 77' meson masses at least marginal 
agreement within errors. However, one should keep in mind the potentially large systematics af- 
fecting the r\' mass determination. 

5.3 Comparing OS and Unitary Approach 

It is also interesting to compare OS and unitary approach on a correlator level. For this we 



computed the correlator matrix Eq. (3.1 ) for D45.32.sc on the same set of configurations with iden- 
tical number of stochastic samples, namely N s = 24. For both cases the GEVP is solved and the 
eigenvalues with errors are exteacted. We then compare in the left panel of Figure [| the relative 
errors of the eigenvalues corresponding to the 77 state X-q by plotting the ratio AA° s /AA,j as a func- 
tion of t/a. This ratio is slightly larger or close to 1. However, the same ratio for the effective 
masses shown in the right panel of Figure || is from t/a = 5 on equal or smaller than 1/2. 

This apparent contradiction is explained by a correlation between A° s (? /a) and A° s (f/a + 1) 
very close to one for all t/a in the OS case. For the unitary case the same correlation coefficient is 
smaller or equal to 0.8. For this reason the tj meson mass can be determined with slightly better 
accuracy from the OS analysis than from the unitary one. The reason is likely to be the variance 
reduction trick for OS strange quarks, as the 77 has a larger strange than light contribution. 



6. Summary and Outlook 

In this proceeding contribution we presented an update of our investigation of r\ and 77 ' meson 
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properties using Nf = 2+1 + 1 flavour of Wilson twisted mass fermions, first published in Ref. [Q]. 
We added one additional ensemble at the finest available lattice spacing. This allows us to estimate 
systematics from the continuum extrapolation with higher confidence. The final result for the 77 
meson mass at the physical point stays virtually unchanged compared to Ref. [^] 

M n =558(14) stat (45) sys MeV. 

In addition we presented an exploratory investigation with a mixed action, where valence strange 
and charm quarks are regularised as so called Osterwalder-Seiler fermions. For one ensemble 
D45.32sc we matched valence and unitary actions using M^. The obtained OS bare strange quark 
mass is in good agreement to the one obtained using kaons for matching. 

It turns out that in the OS approach the 77 can be determined with slightly higher statistical 
accuracy compared to the unitary case, when the number of inversions is approximately matched. 
However, the 77 ' is still noisy and we cannot determine it with higher accuracy than in the unitary 
case. The OS 77 mass shows a smaller dependency on the OS kaon mass than the unitary 17 mass 
on the unitary kaon mass, which we attribute to significant sea quark contributions to the 77 meson. 

In a next step we plan to improve the 77' meson mass determination. This could be reached 
by using the point-to-point method described in Ref. []|] and/or by increasing our operator ba- 
sis. Moreover, we will study the continuum extrapolation of the OS 77 meson mass for different 
matching conditions. 
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